SECOND-ORDER DIFFERENTIABILITY FOR SOLUTIONS OF ELLIPTIC 

EQUATIONS IN THE PLANE 



VLADIMIR MAZ'YA AND ROBERT MCOWEN 



Abstract. For a second-order elliptic equation of nondivergence form in the plane, we in- 
vestigate conditions on the coefficients which imply that all strong solutions have first-order 
derivatives that are Lipschitz continuous or diffcrcntiable at a given point. We assume the 
coefficients have modulus of continuity satisfying the square-Dini condition, and obtain addi- 
tional conditions associated with a dynamical system that is derived from the coefficients of 
the elliptic equation. Our results extend those of previous authors who assume the modulus of 
continuity satisfies the Dini condition. 

Keywords. Differentiability, strong solution, elliptic equation, nondivergence form, modulus of 
continuity, Dini condition, square-Dini condition, dynamical system, asymptotically constant, 
uniformly stable. 

Let us consider an elliptic equation in nondivergence form 

(1) a(x, y) u xx + b(x, y) u xy + c(x, y) u yy = in Q, 

where is an open subset of R 2 . Suppose that u £ W 2 ' 2 {VL) is a strong solution of ([1]); we want 
to know how regular u is in Vt. This, of course, depends upon the smoothness of the coefficient 
functions a, b, and c. If we only assume the coefficients are bounded, then u has first-order 
derivatives that are Holder continuous in f2, i.e. u £ C 1,a (f2) where a £ (0,1) depends on the 
coefficient bounds and the ellipticity constant (cf. 0); if the coefficients are continuous in ft, then 
u e C 1,a (f2) for all a e (0, 1) (cf. 1 ]). On the other hand, if the coefficients are Holder continuous 
in f2, or more generally if the coefficients are Dini-continuous in f2, then it is well-known that 
u 6 C 2 (0) (cf. [6 ). In this paper, we want to find conditions on the coefficients, weaker than 
Dini continuity, under which u will be second-order differentiable. 

Let us assume that f2 contains the origin = (0,0) and focus on the differentiability at 0. 
Using ellipticity and a change of independent variables, we may arrange that a(0) = 1 = c(0) 
and b(0) — 0. Consequently, we assume that the coefficients a,b,c satisfy 

(2) sup (|o(x) - 1| + |6(x)| + |c(x) - 1|) < w(r) as r -> 0, 

|x|=r 

where x = (x, y) and the modulus of continuity w is a continuous, nondecreasing function for 
< r < 1 satisfying to(0) — 0. The coefficients being Dini continuous means that ([2j holds with 
cj(r) satisfying the Dini condition J Q w(r)r _1 dr < oo; Holder continuity, of course, corresponds 
to the special case u(r) — Cr a where a € (0, 1) and C is a positive constant. But we shall assume 
w(r) satisfies the more general square-Dini condition: 



(3) 



f 1 2/ \ dr 

/ lj (r) — < oo. 
Jo r 
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Given a solution u € W 2 ' 2 (il) of let us introduce the vector U = (Ui,U 2 ) = (u x ,u y ). 
Using tin = (Ui) X ) u XJ/ = and w ra = (E72)», we can write jT} as 

0(1, y) (C/i), + 6(1, y) (C/i) y + (c(x,y) - l)(U 2 ) y + (U 2 ) v = 0. 

If we differentiate this with respect to x and use (U 2 ) yx — u yyx — u xyy — (Ui) yy (where third- 
order derivatives are interpreted weakly), we obtain 

(4) (a(x, y) (JJ\)x)x + (b(x, y) (^) w ) a + ((c(x, y) - l)(U 2 ) y ) x + {Ui) m = 0. 

Now we perform a similar calculation using u xy = (U 2 ) x instead of (U±) y and differentiating with 
respect to y instead of x to obtain 

(5) {U 2 ) xx + ((a(x, y) - l)(U{) x )y + (b(x, y) (U 2 ) x ) y + (c(x, y) (U 2 ) y ) y = 0. 
Putting ((4]) and ([5]) together as a second-order system, we obtain 

,» ((; ?H + ((<-> MA$ c °»A(i 

Now may look more complicated than (JT|), but at least it is in divergence form: 

(7) (A n U x ) x + (A 21 U x ) y + (A 12 U y ) x + {A 22 Uy)y = 0, 

where the Ay are (2 x 2)-matrices and U £ W 1,2 (Vt : R 2 ) is a weak solution. Moreover, the matrices 
Aij are perturbations of dijl, where / is the 2x2 identity matrix, in that 

(8) sup |Ajj(x) — Sijl\ < ui(r) as r — > 0. 

|x|— r 

In this way, (O is reminiscent of our work [3 which considered the first-order differentiability of 
weak solutions to an elliptic equation in divergence form. Moreover, the first-order differentiability 
of U corresponds to the second-order differentiability of u, so the conclusions of [3] are just what we 
need here. However, the formulas of [3] pertain to equations and not systems (whose coefficients 
are matrices so that products do not commute); consequently, they cannot be used directly in 
the present situation. Nevertheless, we can apply the methods of [3] to 0. 

The method of [3] suggests that we find a first-order dynamical system on < t < 00 whose 
stability properties as t — > 00 control the differentiability of the solutions of ((?]). To derive the 
dynamical system, we first write x = r9 where r = |x| and 9 = {9\,9 2 ) — (cos (j), sin <f>) for 
< 4> < 2tt. Then we write 

(9) U{x,y) = U (r)+V 1 (r)x + V 2 (r)y + W(x,y), 
where Uq, V%, and V 2 are given by mean integrals 

U (r)=P U{r0)d4, 

(10) J 31 

Vi(r) = f U{re)6xd(t>, V 2 (r) = P U(rd)9 2 d</>, 
J s 1 J s 1 

and W(x, y) has zero spherical mean and first spherical moments: 



(11) / W(r9)d<j) = = Y W{r9)9 l d(j) = y W{r9)9 2 d<f). 

J s 1 J s 1 J s 1 

Similar to |3], we shall show that the 4- vector function V(r) = (Vi(r), V 2 (r)) satisfies a dynamical 
system that depends on W, and W satisfies a PDE that depends upon V. Ultimately, we shall 
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show that the behavior of Uq, V, and W are all controlled by the asymptotic behavior of solutions 
to the following first-order system: 

(12a) -jj- + R<f> = on 0< t < cx), 

where r — e~* and R(e~ t ) is the (4 x 4)-matrix function defined on < t < oo by 

/3i(r) _0 h(r) ci(r)\ 

(12b) R(r):= !fj ^ ° ffj 

a 2 (r) 02 (r) c 2 (r) 

\-ai(r) -6i(r) -ci(r)/ 
with coefficients given by certain second spherical moments of the original coefficients: 

3i(r):=/ a(r6)(9l-6l)d<j>, a 2 {r) := -2 p a(r0)0i02 #, 

(12c) &i(r):=/ b{r6){0l-ei)d<l>, b 2 (r):=-2p o(r0)M 2 #, 

J s 1 J s 1 

ci(r):=/ c(r#)(0 2 2 - 2 ) #, c 2 {r) := -2 f c(r0)M2#. 
J s 1 J s 1 

As in [3], the first-order regularity of solutions of (J7J is determined by the stability properties 
of (|12p . In particular, we say that (|12l) is uniformly stable as f — > oo if for every s > there exists 
a, 5 = 6(e) > such that any solution (j> of (|12ap satisfying |^(ti)| < S for some t\ > satisfies 
1 0(f) | < £ for all t > t\. We shall show that the following holds: 

Theorem 1. If flty) is uniformly stable, then every weak solution U £ W 1,2 (H,,M. 2 ) of Ffy is 
Lipschitz continuous at x = 0. 

Another important stability condition is that a solution of (|12p be asymptotically constant, i.e. 
that (j){t) — > as t — > oo. As discussed in [3], this is actually independent of uniform stability, 
so we need to assume both conditions to conclude the differentiability of weak solutions. We shall 
show that the following holds: 

Theorem 2. If is uniformly stable and every solution is asymptotically constant, then every 
weak solution U € T4 /1 ' 2 (fi,R 2 ) of ^ is differentiate at x = 0. 

Recalling the derivation of (0 from ([TJ, these results yield the following: 

Theorem 3. If US\) is uniformly stable, then every strong solution u £ W 2,2 (fl) of {Ip has 
first-order derivatives that are Lipschitz continuous at 0. //, in addition, every solution of 
is asymptotically constant, then u is second- order differentiable at 0. 

We can obtain analytic conditions on the matrix function R that imply the desired asymptotic 
properties of (fT2|) . The simplest condition is 

(13) r _1 J?(r) £ L^e) for some e > 0, 

which guarantees that (fT2|) is both uniformly stable and asymptotically constant; cf. [2]. 

Corollary 1. If R as in i!2b\) satisfies i!3\) . then every strong solution u £ W 2 ' 2 (il) of (Qp is 
second-order differentiable at 0. 

Analytic conditions weaker than (|13j) can also be obtained. For example, if we introduce the 
symmetric matrix S — —(R+ R t )/2 and let fx(S) denote the largest eigenvalue of S, then it is 
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shown in [3] that 

(14) / p~ l n(S(p)) dp<K for all e > r 2 > n > 



implies that (fl2|) is uniformly stable; as a consequence, (fT4|) guarantees that every strong solution 
u £ W 2,2 (£l) of |lj has first-order derivatives that are Lipschitz continuous at x = 0. In addition, 
it is shown in [3] that 

(15) r^Rir) f p^ 1 R(p) dp € L 1 (0, e) 

Jo 

implies that (|12[) is uniformly stable and asymptotically constant; as a consequence, (JTSJ) guar- 
antees that every strong solution u G W 2 ' 2 (£l) of U) is second-order differentiable at x = 0. 

One may also consider special cases to better understand the significance of the role of the 
dynamical system (|12j) in determining the regularity of strong solutions of (|T|). In particular, let 
us assume that the coefficients b and c in ([l} satisfy 

(16) h(r) =c % {r) = fori = 1,2; 

this occurs, for example, when b and c are constant, or more generally if they depend only on r: 
b = b(r) and c = c(r). In the case ([TB]) . we see that (|12a|) decouples into three scalar equations 

d<f> _ d(f> _ dtp _ 

— + a\ 4> = 0, — - ai 4> = 0, — +a 2 = O. 

dt dt dt 

But these are all of the form <fr' + p(t)<p — which can be solved using the integrating factor 
exp[J p(r)dT]. We conclude that the three scalar equations will be uniformly stable provided 



(17) 



J <xi(t) dr < K\ and J a-zij) dr > —K-2 for t > s sufficiently large. 



Moreover, the three scalar equations will be asymptotically constant provided 

/oo 
3i(t) dr converges to a finite real number, and 

/ C2(t) dr converges to an extended real number > — oo. 
Jt 

Thus, when the coefficients 6, c satisfy (fTB)) and the coefficient a satisfies (|17l) and (TT5)) . then every 
strong solution of (fTJ) will be second-order differentiable at 0. 

1. Derivation of the Dynamical System 
A weak solution U of ^ satisfies 

(19) / Aij diU d 3 ri dxdy = 0, 

Jo. 

for all r\ € Cq°(CI), where we have used the Einstein summation convention of summing over 
repeated indices. To obtain the dynamical system (|12[) . we begin by considering (|19[) with different 
choices of test functions rj. 

Taking rj to be a radial function rj(r) and using ([9]), we obtain the following first-order ODE: 

(20) A(r) U' + rBi{r) V( + rB 2 {r) V{ + r x (r) Vi + T 2 (r) V 2 + A[W](r) = 0, 
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where the (2 x 2)-matrices A, Bk, and Tk arc 



A(r) 



Aij(r6 



(21) 



Si 



I 2 + 0(u(r)) as r -» 
= 0(w(r)) as r — > 



r*(r) 

and the 2- vector A is 
(22) 

Using Lemma 1 in [3J, we can show that 



A ik (r9) Qi d(j> = 0{u(r)) as r -4 0, 



A[W](r) = / (A l i6 i W x + A %2 B % W y ) , 
s 1 



|A[W](r)| < w(r)/ |W| #. 

s 1 



(23) 



Note that, although we are thinking of ([20]) as an ODE, the coefficients are matrices and so it is 
really a system of two equations. 

Taking r\ = r)(r)x and then r\ = Tj(r)y in (|19p . we obtain two second-order ODEs which we can 
put together as a second-order system: 



(24) 



2 /YBi C/A AA U ,4 12 \ /rF/\ /B n B 12 \ fvA (Pi[VW}\ 
B 2 U' ) + \A 21 A 22 ) \rVi) + \B 21 B 22 ) \V 2 ) + \P 2 [VW]J 



r 2 ^ 



621 £22 J KrViJ \C 21 C 



V 2 



Q 2 [VW\ 



0^ 



where the Bj are defined above, but the other (2 x 2)-matrices are 



(25) 



A k e{r) 
B ke (r) 
B u (r) 



Ai^rd) 9i6j6 k e t d<j) = )- S M h + 0{w{r)) as r -> 0, 
4«(r0) 0;0 fc d<?(> = - S u h + 0(uj{r)) as r -> 0, 
A fcl (r<9) 0j0 £ = - <5 fc £ P> + 0(w(r)) as r ->• 0, 
A fc£ (r0) dcf) = 8 u h + 0(uj(r)) as r 0, 



A fei (r0) 0j d</> = 0(w(r)) as r -4 0, 



and the 2-vectors Pi, P 2 , Qi, Q 2 are given by 



(26) flfetVW](r)=> Aiji^OiOk—dcl) and Q fe [V^](r) = / A fel (r6 
J s 1 "^j J , 



dW 



As with ([23]) . we can show 
(27) |A[W](r)|, |Q fc [VW](r)| < w(r)/ \VW\d<f>. 



We want to use (|2T)1) to eliminate C/q from (f2"4")l and then identify the leading-order terms. Since 
A(r) = I n + O(oi(r)) as r — > 0, -4,(r) is invertible for small r and we can write 

BiU'o = -B.A- 1 {B x rV[ + B 2 rV 2 ' + T1V1 + T 2 V 2 + A[VW}) . 
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But the coefficients of rV- and Vj in this expression are "lower-order" , i.e. 

BiA~ x B j} BiA^Tj = 0{uj 2 (r)) as r -> 0. 

So when we plug this into (f24f , it does not affect the leading order terms in rVJ and Vj . Similarly 
for replacing TiU' in 

Now let us make the substitution r = e - *, so that rd/dr = —d/dt. Next, let us introduce 

e{t) = u(e-% 
and then write (|2~4")1 (after the elimination of Uq) as 

e~ 2t (-AV t + BV + P[VW] + 0(e 2 (i))^ + e - 2t (-B V t + C V + Q[VW] + 0(e 2 (t))) =0, 

where V — (V\, Vz), P = (Pi,P2), and Q = (Qi,Q2) are 4-vectors and A, B, B, and C are the 
(4 x 4)-matrices 

A = A 2 )- b = (b" 5"). § = (f" !")• c =(e" c" 

\Si2i -A22J \&21 £>22/ \021 O22) \ L 21 ^22 

and we have used 0(e 2 (t)) to represent terms depending linearly on Vt, V, or A[VW], but with 
coefficients that are 0(e 2 (t)) as t —> 00. We can remove the factor e~ 2t to obtain 



(28) 



-AV t + BV + P[VW] +0{e 2 {t)) +(2A-B)V t 

i t 

+(C - 2B)V - 2P[WW] + Q[VW] + 0{e 2 {t)) = 0. 

However, this is still a second-order system, and we want to avoid differentiating the coefficient 
matrices, so let us convert it to a first order system by replacing the vector in the brackets in (|28[) 
by a new 4- vector 

(29) U = -AVt+BV + P[\7W]+0{s 2 (t)). 

We now have a first-order system in the 8- vector (V, U): 

V t - A X B V + A' 1 U = A~ l P[VW] + 0{e 2 ) 
U t + (C- BA X B)V + (BA- 1 - 21)U = -Q[VW] + 0(e 2 ). 
where I is the (4 x 4) identity matrix. The coefficients of V and U behave as follows: 

-A 1 B~-I, A" 1 ~2I. C BA *B ~ -I, BA 1 21 - I. 

2 ' 

where ~ means differs by a term that is 0(e) as t —¥ 00. Consequently, let us rewrite the first-order 
system as 

where the (8 x 8)-matrix-valued function M(i) is of the form 

M(t) = Moo + Si(t) + S 3 (t), 

with a constant matrix 

The variable coefficient matrices Si and S2 satisfy 

= (c -BA-^- |I BA-^- 2 l) =0 ^ aS< ^°° 
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and S2 = 0(e 2 (t)) as t —> 00. The right-hand side of ([30]) satisfies 

\Fi(t,VW)\ <e(t)f \VW\dcf>. 
J s 1 

In order to analyze ([50]) . as in [3] we introduce a change of variables 
(31) 

where the matrix 



u) ~" \ > r 



21 21 
1 1 

diagonalizes M^, i.e. J _1 M 00 J = diag(0, 0, 0, 0, —2, —2, —2, —2). We find that (</>,V0 satisfies a 
dynamical system of the form 



where 



with 



R(i) 



i?i(i) i?2(*) 



iJx(t) « ^A" 1 - ^A _1 B + C - BA _1 B + ^BA" 1 - I, 
where « means differs by a term that is 0(e 2 (t)) as t — ► 00. The right-hand side of satisfies 



(33) |G(f, W)l < e(f) / |W| ^ 

Estimates on that we shall discuss in the next section together with the stability theory 
presented in Section 2 of [3J show that the stability of (|52"j) is determined by that of 

(34) ^ +jRl (^ = , 

so we need to determine the asymptotic behavior of R\. But to do this, let us write 
A = i(I + A ), B = i(I + B„), B = i(I + B ), C = I + C , 

where |A |, |B |, |B |, |C | = 0{e[t)) as t 00. Also note that 

A" 1 w2(I-Ao). 
Using these we can simplify R\ to obtain 

i?i ~ Co — Bo = C — 2B, 
and after a careful calculation we obtain the formula given in (|12bl) . 
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2. Proofs of Theorems 1 and 2 

Since we are only interested in the behavior of our weak solution near 0, we may assume 
ft = B £ (0) with e > chosen small enough to make 

(35) / r _1 u;(r) dr < S and u(e) < 5, 



with S > as small as we like. In fact, for any p G (1, oo) we can choose 5 = 8(p) > in (|35|) small 
enough that the small oscillation condition on the coefficients ([2]) ensures that Vt/ € L^ oc (ft); cf. 
Corollary 6.2 in [4]. Henceforth, we pick p > 2 and choose e small enough that VZ7 € L£ OC (ft). 
But by rescaling the independent variables, we may arrange e > 1, so we may assume that our 
weak solution U of {7J satisfies 

(36) V£7 G L p (ft) where p > 2 and ft = £i(0). 

In particular, by Sobolev's inequality we know that U is continuous in ft. 

For our analysis, it is useful to consider (0 on all of M 2 , so we extend the matrices Ay to all 
of M 2 by 

Aij = Sijl for |x| > 1. 

We also extend our modulus of continuity w to (0, oo) by w(l) for r > 1. It will also be useful to 
introduce the L p -mean of a function over the annulus A r = {x : r < |x| < 2r}: 

i/p 



M P (f,r) = \f(x)\?dx 



To control growth of the first derivatives of functions, we introduce 

M hp (f, r) - rM p (V/, r) + M p (f, r). 

Let us introduce a smooth cut-off function x( r ) t na t is 1 for < r < 1/4 and for r > 1/2. 
We find that x( r )U(x,y) satisfies 

{Au(xU) x ) x + (A 21 ( x U) x ) y + (A 12 (xU) y ) x + {A 22 ( X U) y )y = F + {F r ) x + [F 2 ) y 

where 

F = A llX 'e x U x + A 2lX '8 2 U x + A 12X 'e x Uy + A 22 x'9 2 U y , 
Fx = + A 12 9 2 )U and F 2 = x'iMiOiU + A 22 9 2 )U. 

Using p^l> with r] = X , we see that 

(37) / F dxdy = 0. 



Since we are interested in the behavior near x = = y where U and xU agree, we can simply 
assume that U is supported in r < 1/2 and satisfies 

(38) d i (A ij d j U) = F +d i (F i ), 

where F ,F 1 ,F 2 G L P (R 2 ) are supported in 1/4 < r < 1/2 and F satisfies J37]). Of course, we 
now must replace (fT9| by 



(39) / A^ djUdir] dxdy = t (FAv- F )dxdy, 

JR 2 JR 2 

for all 77 G C^°(ft). 

At this point we observe that (|38)) with p7|) for the vector function U is identical with (51ab) 
in [3] for the scalar function u. This means that we can repeat the analysis of [3] to connect the 
stability of the dynamical system (|12bj) with the regularity of our weak solution. We do not want 
to repeat all of the details here, but let us give an outline of the argument. 
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To begin with, we recall the decomposition U = Uo + Vxx+V 2 y+W in (|9|). We have shown that 
V satisfies a dynamical system (f5U|) that depends on VW, so we need to know VW is sufficiently 
well-behaved in order to obtain estimates for V. This is done by showing that W satisfies a PDE 
that depends on V. To derive the PDE for W, we introduce 

(40) Qij = Aij - SijI, 
which satisfies 

|fiy|<u;(r) for < r < 1 and fly = for r > 1. 
We also introduce for / e L] oc (M 2 \{0}) 

(41) /(r0) x = f(r9) - Pf(r8), 

where P is the projection of / onto the functions on S 1 spanned by 1, 6\, 9 2 : 

Pf(r9) = co(r) + c x (r)0i + c 2 (r)9 2 , where 

Co(r)=/ /(r0)# and Ci {r)=2p 9,f{r9)d<t>. 
J s 1 J s 1 

Notice that 

P[A(U a + Vix + V 2 y)] = A(U + V x x + V 2 y) and P[AW] = 0, 
so W satisfies the following perturbation of Laplace's equation on R 2 : 

(42) AW + [O^djUo)] 1 - + [d l (a iJ d ] {V k x k ))] 1 - + [diiVtijdjW)} 1 - = [F + d^F,)]^ . 

Now we simultaneously consider the dynamical system (|5P|) for V and the PDE (l4"2l for W. 
The analysis in [3] shows the assumptions that U £ W 1,2 (£l) and that (|T2|) is uniformly stable 
together imply that V satisfies 

(43) sup (\V(r)\+r\V'(r)\) < C 

0<r<l 

and that W satisfies 

(44) M hp (W,r) < Cuj{r)r for < r < 1. 

(In both (j43]l and (j44|) the constants C depend upon the W /1 ' 2 -norm of U, but not on r.) Since 
p > 2, we can use Sobolev embedding to conclude that |W(a;, y)^ 1 < Cui(r) for < r < 1, i.e. 
|W(x, — > as r — |x| — > 0. This shows that W is differentiable at 0. 

To estimate Uq, we use (p0|) and the estimates that we have obtained on Vi and VW to conclude 



|f/ o (r)-f/ o (0)| 



<Cuj(r) / (pIF'^I + I^I + IVWI) dp 
Jo 

< Cw(r) r. 

But this implies that f/o is differentiable at x = and ?7q(0) = 0. 



We have now shown the assumption that (fl2|) is uniformly stable is sufficient to show that our 
weak solution U € W 1,2 (f2) satisfies 

\U(x,y) - (7(0)1 < \U (r) - U o (0)\ + |Vi(r)| • |x| + |V 2 (r)| • |y| + |W(x,y)| < Cr, 

But this shows that U is Lipschitz continuous at x = 0, completing the proof of Theorem 1. 

For Theorem 2, we add the assumption that every solution of (|12[) is asymptotically constant. 
The dynamical systems analysis of [3] applied to (13"2l then shows that <j>(t) — > tfioo and i/j(t) — > 
as t — > oo. However, we can use (|3ip to express <j), ip in terms of V and Vt: 



0(e(t)). 
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Hence the conclusion ip — » implies Vt — ► as t -> oo, in other words 
(45) lim rV'(r) = 0. 

But (|45l implies that V\{r)x + Vi(r)y is differentiable at 0. Since we have already shown that Uo 
and W(x, y) are differentiable at 0, we obtain the conclusion of Theorem 2. 
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